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I. INTRODUCTION AND OVERVIEW 


System identification technology has been used successfully 
for many vehicles. Because of their large number of degrees of 
freedom and complex aerodynamic interactions, the rotorcraft have 
always presented a special challenge to system identification 
methods. A completely integrated methodology has been developed 
under this NASA contract to solve this difficult problem. This 
methodology has also been translated into a user oriented series of 
computer programs. This volume provides basic guidelines for effi- 
cient and effective use of one of these computer programs. 

Figure 1.1 shows a schematic flowchart of the overall data pro- 
cessing technique for rotorcraft. The first step in this procedure 
is state estimation and instrument calibration. This is implemented 
by the computer program DEKFIS (for Discrete Extended Kalman Filter 
and Smoother) which implements an extended Kalman filter/smoother 
using the Friedland-Duffy formulation. Instrument biases and scale 
factors are estimated at this stage together with any state which 
is not measured directly. The second step involves estimation of 
the mathematical model of various forces, moments and interchanges. 
This is implemented in OSR (Optimal ^ubset Regression) computer pro- 
gram which uses a regression technique. Accurate estimates of 
parameters are obtained in the final step. One of two computer pro- 
grams is used for this purpose. SCIDNT implements the maximum like- 
lihood method for linear systems and NLSCIDNT extends the method to 
nonlinear rotorcraft models. 

Accuracy of parameter estimates may be improved by using flight 
test inputs based on the input design program, INDES. 

This user’s manual describes the INDES computer program. The 
details of the theory and the particular implementation used are 
given in the final report.* 


Hall, W.S., Gupta, N.K., Hansen, R. and Bohn, J., "State Estimation 
and Parameter Identification for Rotorcraft," final report on con- 
tract NASI-14549, May 1978. 
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II. PROGRAM DESCRIPTION 
2.1 PROGRAM STRUCTURE 

The overall logic of the INDES program is presented in 
Fig. 2.1. It is intended to give the user a general overview of 
the program and, therefore, it shows only the major routines with 
a brief description of their purpose. 

The subroutine calling structure is given in Fig. 2.2. The 
user may find this figure useful in further understanding the pro- 
gram's flow and in constructing an overlay structure if one is 
needed. The program was developed on a CDC 7600. 

The functions of the various subprograms are sketched below: 

MAIN is the main routine. It performs no claculations itself, 
but rather calls three routines --INPUT, INDES, and 0UTPUT- -which 
accomplish the computational and input/output tasks. MAIN's most 
important function is the dimensioning of the major arrays and set- 
ting up most of the labeled commons in the program. 

BLKDAT is a block data routine. Its major function is to 
initialize arrays which contain labels for the measurement, control, 
and expansion variables for input/output purposes. 

INPUT handles all of the program's input, checks for errors in 
the inputs, and initializes various flags within the program. It 
also prints this information as a record for the user. 

FLAG IT is called by Subroutine INPUT to set a pointer arrays 

to the state, measurement, and expansion variables which are to be 
used in the run. 

SYMWl is called by INDES to compute the eigenvalues and 
eigenvectors of the (symmetric) information matrix. It is optimized 
for use on symmetric matrices. The eigensystem is calculated by 
the QR algorithm after the matrix has been put in tridiagonal form 
by a Householder transformation. 
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Figure 2. 2 Subroutine Calling Structure 




















UPDATE ; (1) propagates the state estimates (x) and state sen- 

sitivities (3x/39) with respect to the parameter estimates and com- 
putes the measurement estimates (y) and their sensitivities (3y/36) 
by calls to appropriate subroutines; and (2) stores the measurement 
estimate time histories for use by 0UTPUT. 

INTGR8 is called by UPDATE to oversee the integration of the 
first-order differential equations for x and 3x/30 forward one 
sample time increment. This is accomplished by a variable-order, 
variable-stepsize, predictor-corrector algorithm (Adams -Krogh Method) . 
The algorithm chooses its own order and stepsize to most efficiently 

satisfy relative and absolute error bounds on x and ax/ae set by 
the user. The resulting stepsize may be more as well as less than the 
sample time increment. 

STEP is called by INTGR8 to perform the actual integration of 
the equations. In general, the integration stepsize is not a 
rational fraction of the sample time interval. 

INTRP is called by INTGR8 to interpolate the values of x, 

3x/3§ and their time derivatives at the data sampling times. 

DXDTNL is called usually by STEP but initially by INTGR8 to 
obtain the time derivatives of x and also 3x/3§^, where 
is an identified parameter. 

YAHAT is called by UPDATE to evaluate y and 3y/36 at the 
data sample times given x and 3x/36 at those times . 

STATE has five ENTRY points. A call to STATE initializes var- 
iables used elsewhere in the subroutine. STATIC sets the state 
initial condition estimates, x(0) . STAICI sets the state sensi- 
tivity initial condition estimates, 3x(0)/38. A call to XD0T com- 
putes the time derivate of the state estimates, dx/dt. A call to 
XTHD0T computes the time derivative of the state sensitivities with 
respect to a specified parameter, d(3x/30^)/dt. Therefore, the 
vehicle equations of motion appear in this subroutine. 
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MEAS has three ENTRY points. A call to MEAS evaluates the 
measurement noise covariance matrix, R. HC0MP computes the meas- 
urement estimates, y. HTHCMP computes the measurement sensitivi- 
ties with respect to a specified parameter, 3y/30^. Therefore, 
the measurement instrument equations appear in this routine. 

C0EF evaluates the aerodynamic coefficients required by the 
vehicle equations of motion in subroutine STATE. 

DC0EF evaluates the partial derivatives of the aerodynamic 
coefficients with respect to each of the identified parameters. 

C0NTRL finds the values of the control inputs at any specified 
time by linear interpolation between the nearest sample time points. 

0UTPUT performs one major task: it uses PLTLCM to plot the 

final control input time histories and to plot the resulting simu- 
lation. 

In addition to these major subroutines, there are 15 utility 
routines of which use is made by several of the above routines. 

ADD finds the sum of two matrices. 

EQUATE sets one matrix equal to another. 

INV finds the inverse and determinant of a matrix. 

INV2 finds the inverse and determinant of a matrix. 

MULT finds the product of two matrices. 

PRNT writes a matrix on the printer. 

TRCE computes the trace of a matrix. 

RDTITL reads a card containing the run title. 

TRANP tranposes a matrix. 

SUBT finds the difference of two matrices. 

ZER0 sets a matrix equal to the zero matrix. 

ASPERR produces a subroutine trace -back in the output when an 

error in computation is detected. 
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1 

I LNCNT prints the run's title at the top o£ each page. 

PLTLCM produces plots on the printer. 

J 2.2 INPUT DESIGN ALGORITHM 



The following steps detail the computational procedures indi- 
cated in the various boxes of Fig. 2.1. 


2.2.1 The nominal input i “ 1,N*NU is given (input) 

N*NU 



• k is input quantity 

• llnew, ^ (N*NU)*1 vector 


where is a perturbation of h^^ — scaled by /k, typically, 

k=.01. 

2.2.2 h is composed of NU vectors of length N stacked 

new 

together. This can be considered as NU separate vectors 
of length N with i'=l,NU indicating the vector and 
ilocai“i>^ indicating the component. 

2.2.3 y^ for each £=1,NM and for each ^*1,NM and 

k=l,NP are simulated by the fipDATE routine in NLS. 
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In NLSC 


YACt) 


I(t) 

iz(t) 

36j^ 


3l(t) 

3®NP 


Output response 


Output sensitivites 


to 0. 


'NP 


2.2.4 The time histories are to be averaged over each interval so 
that YACt) can be expressed as N "discrete” values for each 
time history 


2.2.5 The nominal output response and sensitivity function time 
histories are stored in the d array. 


2.2.6 To peturb at the ith position by 6p (which is an 

input quantity) 


►NLOCAL = 0 
-DO i ' = 1,NU 


DO i 


LOCAL 


= 1,N 


ip = (i'-l)*N + 

^new^^p^ “ ^new 
pert 

^new^^p^ “ ^new 
pert 


ONTINUE 

CONTINUE 


^LOCAL 

(ip) (for ip 


^ i3 
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N NM 

2.2.7 Q-Ci,,!?) = 2 2 pCinJ.l.i)- pCi2 

^ j«l £=1 J- ^ ^ 

+ i Q-(((L,-l)/N+n 

Note: Cl) the ’'1" index in p is for the output time history 

which is assumed to be stored before the sensitiv- 
ity function time histories. 


(2) 6i,i 


1^2 = (0 if ii f* i 
U if ii * i 


2.2.8 To calculate M, 


N*NU 


Cl) ReCj»h,5,) = Z h Cm) pCm,j,k,£) for all k>l 
^ m=l 

• N NM 

C2) M^Cki-l.k^-l) = Z Z CRJj.k, £) • RJj,k~ £)/l^Ci)) 
s 1 2 : j=l £=1 ^ s 


2.2.9 To calculate 17 Crequires inverting M,,) 
For Case 1 


i7 = CM'^)^ W 

1 CW is input) 
J3 =• T^{W M3I) 


For Case 2 


= M. 


-1 

1-1 


J3 = im;*i 


2.2.10 Calculate V’ 

NP+1 NP+1 

vein I 2 ) = 2 ^ I 

k,=2 k.,=2 


N NM 

WCk,-l,k2-l) -2 2 

^ ^ j=l il=l 


Pdi^ j ,k^^£)pCi2^J ) 

RC£) j 


2.2.11 Want to solve the generalized eigenvalue problem: 

V'h = XQ’h 

for the largest eigenvalue X and the corresponding eigenvector h. 
V and Q' are both symmetric. 
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The IMSL routine EIGZF is used which obtains all the eigen- 
values for the general problem. For a thorough explanation, see 
the IMSL documentation, then sort through these for the largest 
eigenvalue and the corresponding eigenvector. 


2.2.12 To calculate ^stop 


For Case 1 - Tr {W 

For Case, 2 X^^^p - NP » 


2.2.13 /6 ® - /ay + /ay^-a+l 

(where is input) 

h^ = /a hg + /F h 

2.2.14 Y * hg’^Q'h 


2.2.15 To calculate M 

N*NU 

(1) R Cj,k,Z) = Z h(m)p (m, j ,k, a) for all k>l 

m»l 

N NM 

(2) M (k,-l,k2-i) = Z Z (R(j,k, i) • RCj.k. iO/F(il)) 

^ j=l £=1 

Note that M is a NP*NP matrix. 


2.2.16 To calculate J (must invert M) 

For Case 1: J = Tr {W M’^} 

For Case 2: J = lM"^| 


2.2.17 To calculate /6 


Solve a + 6+2 /a6Y=l 


for 0 </6^ 1 
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i.e. ^ - /a y * /ay^-a+l 

2.2.18 h' = /5T hg + /6 h 

2.2.19 To calculate M' 

N*NU 

(1) R_Cj,k,il) = Z h (m) p(m,j,k,Jl) for all k>l 

® m=i ^ 

(this should be available from 2.2.8) 

N NM 

(2) M*(k,-l,k2-l) » Z Z RJj,k, JO • R(j,k. Jl)/l?(JO 

^ ^ j=l il=l ® 

(3) M' = a Mg + 6M /oB M* + /o^ (M*)”^ 

2.2.20 To calculate J* see 2.2.16 above -- invert M* and use it 
instead. 



Other comments: LIMITJ | 

> are input on card types 6 and 4 
J respectively (see Table 3.1) 


13 




III. PROGRAM INPUT 


3.1 CARD INPUT 

Table 3.1 gives the details of how to set up the input data 
cards. Tables 3.2 through 3.4 give additional information needed 
to generate the input data and are referred to in Table 3.1 when 
the information is needed. The parentheses in Table 3.4 indicate 
parameters of secondary importance which are usually left unused. 

A listing of input data is shown in Figure 3.1. 


IS 



Tabl 
INDUS Input 


Cols. 

Format 

1-80 

15A4 

1-5 

15 

6-10 

15 

11-15 


16-20 

15 

21-25 

15 

1-20 

415 


1-10 F10.0 


Name 

TITLE 

N 

NU 

NS 

NM 

NP 

NC(I). 
1=1 ,NU 

T 


21-30 


FIO.O 


DELTAP 


3.1 


Card Setup 


Description 


Alphanumeric title for problem 


Number of time steps in the designed control 
input sequence 

Number of simultaneous control inputs to be 
designed 

Number of states in the assumed model (see 
card type 15) 

Number of measurements in the assumed model 
(see card type 15) 

Number of parameters in the assumed model 
that are to be identified (i.e. the input is to 
maximize the identif lability of this many 
parameters; see card type 16) 


Vector indicating which of the eight control 
inputs in the model are to be used. The 
value of NC(i) is the appropriate index (see 
Table 3.2-Control Variables) for the ith 
control in the simultaneous design. 


Total time in maneuver - should be compatible 
with NN and DELTA in card type 5. I.e. 
T=(NN-1)*DELTA 


Magnitude of perturbation from nominal input 
to be used in calculating sensitivities to 
control step changes. 













Table 3.1 (Continued) 


Number of 
Type Cards 

Cols. 

Format 

Name 

Description 


31-40 

FIO.O 


Allowable error for computing convergence 
criteria on largest eigenvalue (see Figure 2.1) 


41-50 

FIO.O 


Allowable error for computing convergence 
criteria for J loop (see Figure 2.1) 


51-60 

FIO.O 

EE 

Use 0.05 (other options not Implemented) 

1 

1-5 

15 

NN 

Number of time points In continuous time 
history simulation 


6-15 

FIO.O 

DELTA 

Step size in continuous time history 
simulation 

1 

1-5 

15 

NCASE 

Case type: NCASE =1 J = TR 

NCASE = 2 J = |M-1| 


6-10 

15 

LIMITI 

Max. number of I loop Iterations (see Fig. 2.1) 


11-15 

. 

15 

LIMITJ 

Max. number of J loop Iterations (see Fig. 2.1) 


16-20 

15 

LIMITO 

Use 1 (other options not Implemented) 


21-25 

15 

KITLIM 

Use 3 (other options not Implemented) 

NU- 

Alternately 
spaced with 
an equal 
number of 
cards of 
type 8, I.e. 
7-8- 7-8. . . . 

1-80 

8F10.0 

HN(i). 

i.=l,M 

This Is the nominal nonlinear control Input 
about which an optimal perturbation Input is 
to be designed. This Is typically a constant 
or zero input setting, but can also be a 
maneuvering sequence. The optimal control 
Input Is the sum of this nominal and the 
optimum perturbation control Input that Is 
designed. 

NU-see note 
on card 
type 7 

1-80 

8F10.0 

HS(i), 

1=1 ;m 

This is the initial guess for the perturba- 
tion (from nominal) control Input that Is to 
be optimized 




















Table 3.1 (Continued) 


Card 

Type 

Number of 
Type Cards 

Cols. 

Format 

Name 

Description 

9 

(k*NP), 

1-80 

8F10.0 

i=l,NP 

W weighting matrix for the Identification para- 
meters, 1 row per card (used for Case 1). 
k=[NP/8] truncated +1 

10 

1 (or 2 as 
necessary) 

1-80 

8F10.0 

im^ii 

Diagonal elanents of the noise covariance 
matrix 

11 

1 

(or 2 as 
necessary) 

1-80 

8F10.0 

i=l,NM 

parameter weighting terms on outputs In 
energy constraint expression (see Figure 2.2.7) 

12 

1 

1-20 

8F10.0 

; 

Q 2 ( i ) * 
i=l,NU 

Weighting terms on control Inputs in energy 
constraint expression. 
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1 

1-5 

6-10 

11-15 

16-20 

21-25 

26-30 

15 

15 

15 

15 

15 

15 

lOPRT 

lOPLT 

lOPRTl 

lOPLTl 

IRW 

IRWUNT 

Use 4 (other options not Implemented) 

Use 0 (other options not Implemented) 

Use 0 (other options not Implemented) 

Use 0 (other options not implemented) 

= 1 Writes the control perturbation array 
(P array) to file IRWUNT after 
calculating 

=2 Reads the control perturbation array 

(P array) from file IRWUNT. This avoids 
recalculation 

= 0 Calculates P array but does not write 
to file 

Output file for storage of control perturba- 
tion array (P array). If blank, defaults to 
Unit 2. 






















Table 3.1 (Continued) 



Format 

Name 

Description 

11 

IPL0T 

= 3 if both printer plots and a tape storing 
the time history data are to be made 

=2 if only the tape is to be made 

=1 if only the printer plots are to be made 

=0 if neither 

11 

IPRNT 

Sets level of detail of the diagnostic print- 
out 

14 

INCPR2 

= k if X, 8x/3§, y, 3y/a§, and y are to be 
printed every kth sample point 

14 

INCPR3 

Default value is 50. INCPR3 is ignored if 
IPRNT < 3. 

14 

INCPLT 

= k every kth data point is plotted in the 
printer plots (default =1) 

A2.8X 

LTYPE 

= STATES, if the card lists the input names 
of states which are to be integrated 

= READ, if the card lists the input names of 



the states which will be looked-up from 
input time histories 

= MEASURE, if the card lists the input names 
of the measurement variables 

= EXP VARS, if the card lists the input names 
of the expansion variables 

= *b, if no more cards of this type are to be 
read ("b" indicates a blank). 

Note that only the first two characters are 
actually read; these words must be left 
justified. 










Table 3.1 (Continued) 


Card 

Type 

Number of 
Type Cards 

Cols. 

Format 

Name 

Description 

15 

Contd. 


11-80 

14(A3, 

2X) 

(LL(J). 

J=l,14) 

Input names of state, measurement, or expan- 
sion variables as listed in Table 3.2. Note 
that only the first three characters are 
actually read; these words must be left 
justified. Blank words may appear between 
input names. 

16 


1 

A1 

ECHK 

= blank, if this card contains parameter 
information 






= *, if no more cards of this type are to be 
read 



2-4 

13 

J1 

Parameter index (see Table 3.3) 



5 

A1 

J2 

= *, if this parameter is to be identified 
= blank, if not 



6-11 

A6 

PLABJl 

Label for the parameter for printout purposes 



12-30 

E19.0 

PJl 

Initial value of the parameter 



31-32 

12 

11(1) 

= k, this parameter appears in the polynomial 
expansion of the kth aero coefficient 
(see Table 3.4), 






Ignored if J1 > 400. 



34-50 

17L1 

J=2,..,18 

1 

Exponents of the expansion variables in the 
aero coefficient polynomial term containing 
this parameter 

17 


1 

A1 

JA 

= T, the label which follows on this card is 
for time (the independent variable) 






= Y, the label is for a measurement variable 
= U, the label is for a control variable 






= X, the label is for a state variable (only 
look-up states are affected) 






Table 3.1 (Concluded) 


Card 

Type 

Number of 
Type Cards 

Cols. 

Format 

Name 

Description 

17 

Contd. 

^ 

3-4 

12 

J 

The index number of the variable. If JA = T, 
then J is ignored. 


- 

7-54 

8A6 

1 

NAME 

i 

• ■ 1 

A label of eight 6-character words. 



Table 3 . 2 

State, Control, and Measurement Variables 


SYMBOL 



STATE VARIABLES 


DEFINITION 


longitudinal component of velocity 

lateral component of velocity 

vertical component of velocity 

body roll rate 

body pitch rate 

body yaw rate 

Euler roll angle 

Euler pitch angle 

Euler yaw angle 

coning angular rate 

longitudinal flap angular rate 

lateral flap angular rate 

coning angl e 

longitudinal flap angle 

lateral flap angle 

collective lag angular rate 

longitudinal lag angular rate 

lateral lag angular rate 

rotor speed variation 

collective lag angle 

longitudinal lag angle 

lateral lag angle 

rotor angular position variation 






Table 3.2 (Continued) 


j 


CONTROL VARIABLES 


SYMBOL 

■ DEFINITION 

®0 

collective pitch of blades 

®lc 

lateral cyclic pitch of blades 

®ls 

longitudinal cyclic pitch of blades 

<^TR 

pitch of tail rotor blades 


elevator angle 


aileron angle 

6 

r 

rudder angle 


flaperon angle 


MEASUREMENT VARIABLES 


UNITS 


r<jg 

rad 

rad 

rad 

rad 

rad 

rad 


SYMBOL 

DEFINITION 


longitudinal accelerometer 

^y 

lateral accelerometer 

J 

vertical accelerometer 


roll angular accelerometer 


pitch angular accelerometer 


yaw angular accelerometer 


rol 1 rate gyro . 


pitch rate gyro 

■•ni 

yaw rate gyro 


roll position gyro 


pitch position gyro 

lb 

yaw position gyro 


angle-of-attack vane 

Sfn 

sideslip vane 


pitot tube 


UNITS . 

ft/sec 

ft/sec^ 

ft/sec^ 

rad/sec 

rad/sec 

rad/sec 

rad/ sec 

rad/sec 

rad/ sec 

rad 

rad 

rad 

rad 

rad 

ft/ sec 


2 

2 

2 













Table 3.2 (Concluded) 


MEASUREMENT VARIABLES (Cont'd) 


f^DEX 

SYMBOL 

DEFINITION 

UNITS 

16 

"m 

longitudinal velocity 

ft/sec 

17 


lateral velocity 

ft/sec 

18 


vertical velocity 

ft/sec 

19 

XRm 

rotor longitudinal force 

lbs 

20 

''Rm 

rotor lateral force 

Tbs 

21 


rotor vertical force 

lbs 

22 

'■Rm 

rotor roll moment 

ft/lbs 

23 

% 

rotor pitch moment 

ft/lbs 

24 

'‘Rm 

rotor yaw moment 

ft/lbs 

25 

^Im 

flap angle of blade 1 

rad 

26 

^2m 

flap angle of blade 2 

rad 

27 

®3m 

flap angle of blade 3 

rad 

28 

. ^4m 

flap angle of blade 4 

rad 

29 

^5m 

flap angle of blade 5 

rad 

30 

^6m 

flap angle of blade 6 

rad 

31 ‘ 

^7m 

flap angle of blade 7 

rad 

32 

^Im 

lag angle of blade 1 

rad 

33 

^2m 

lag angle of blade 2 

rad 

34 

'’3m 

lag angle of blade 3 

rad 

35 

^4m 

lag angle of blade 4 . 

rad 

36 

^5m 

lag angle of blade 5 

rad 

37 

^6m 

lag angle of blade 6 

rad 

38 

‘’7m 

lag angle of blade 7 

rad 

39 

♦R>n, 

cosine of rotor azimuth 

N.D. 

40 


sine of rotor azimuth 

N.D. 
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Table 3.3 
Parameters 


INDEX 

SYMBOL 

DESCRIPTION 



These parameters are available for defining 



the polynomial expansions of the aerodynamic 



coefficients 

401 

u(0) 



402 

v(0) 



403 

w{0) 



404 

p(o) 



405 

q(0) 



406 

r(0) 



407 

<J>{0) 



408 

e(o) 



409 

<1^(0) 



410 

e„(o) 



411 

8i (0) 


^ Initial conditions for the state 


Ic' ' 


variables 

412 

8ls«» 



413 

Bo(0) 



414 




415 

fls'O) 


■ . 1 

416 

=0«>' 


• 

417 

he"” 



418 




419 

4ni,(0) 



420 

'o"» 



421 

he"” 



422 

hs"” 



423 

B'I'r(o) 

J 


423-425 


Not used. 

426 

9 

Gravitational acceleration (default - 



9.80665 m/sec^) 

427 

p 

Air density 

428 


Rotor shaft tilt in X-Z plane, (+) forward 







Table 3.3 (Continued) 


INDEX 


SYMBOL 



DESCRIPTION 


Rotor speed 
Rotor radius 
NC/ R 

Lock number = pacR^/I. 

b 

2-0 lift curve slope 
Mass of vehicle 

Longitudinal distance of rotor hub from 
vehicle C. 6. 

Vertical distance of rotor hub from vehicle 
C.G. 


Inertia terms 


Number of blades 


26 




Table 3.3 (Continued) 


INDEX 


SYMBOL 


459 

460 

461 

462-469 


470 

471 

472 

473 

474 

475 


479 

480 

481 

482 

483 

484 


488 

489 

490 

491 

492 

493 


cgx 


cgx 


'cgx 


cgy 


cgy 


'cgy 


cgz 


cgz 


'cgz 


DESCRIPTION 


Inertia terms 


Longitudinal Accelerometer 
Noise covariance 
Bias 

Scale factor error 
X location 
Y location 
Z location 


Lateral Accelerometer 
Noise covariance 
Bias 

Scale factor error 
X location 
Y location 
Z location 


Vertical Accelerometer 
Noise covariance 
Bias 

Scale factor error 
X location 
Y location 
Z location 












Table 3.3 (Continued) 


SYMBOL 


DESCRIPTION 


Roll Angular Accelerometer 


Noise covariance 
Bias 

Scale factor error 
Pitch Angular Acclerometer 


Noise covariance 
Bias 

Scale factor error 
Yaw Angular Accelerometer 


Noise covariance 
Bias 

Scale factor erro 
Roll Rate Gyro 


Noise cavariance 
Bias 

Scale factor error 
Pitch Rate Gyro 


Noise covariance 
Bias 

Scale factor error 
Yaw Rate Gyro 


Noise covariance 


Scale factor error 




i 

1 

Table 3.3 (Continued) 


INDEX 

SYMBOL 

DESCRIPTION 

! 


Roll Position Gyro 

533 

"r 

Noise covariance 

534 


Bias 

535 


Scale facotr error 


Pitch Position Gyro 

539 Rg Noise covariance 

540 bg Bias 

541 kg Scale factor error 


545 

% 

546 ■ 

k. 

il> 

547 

551 

R 

a 

552 

b 

a 

553 

k 

a 

554 


cga 

555 

•^cga 

556 

557 

^cga 

558 

^6 

^6 

^cge 

559 

560 

561 

562 

•^cgB 

563 

^cgB 


Yaw Position Gyro 
Noise covariance 
Bias 

Scale factor error 

Angle-of-Attack Vane 
Noise covariance 
Bias 

Scale factor error 
X location 
Y location 
Z location 

Sideslip Vane 

Noise covariance 
Bias 

Scale factor error 


X location 
Y location 
Z location 








Table 3.3 (Continued) 



INDEX I SYMBOL DESCRIPTION 


Pitot Tube 

Noise covariance 
Bias 

Scale factor error 
X location 
Y location 
Z location 

Reference air density 
Velocity of sound 

Longitudinal Velocity Measurement 


Noise covariance 

575 Bias 

576 Scale factor error 

577 

Lateral Velocity Measurement 


578 Ry Noise covariance 

579 Bias 

560 Scale factor error 

581 

Vertical Velocity Measurement 


562 R^ Noise covariance 

563 b^^ Bias 

564 k^ Scale factor error 

585 

566 Horizontal reference distance 

587 

^ref Vertical reference distance 

589 

Rotor Longitudinal Force Measurement 
500 Rxd Noise covariance 


Scale factor error 
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Table 3.3 (Continued) 




INDEX 

SYMBOL 

DESCRIPTION 



Rotor Lateral Force Measurement 

594 

^yR 

Noise covariance 

595 


Bias 

596 


Scale factor error 

597 

Rotor Vertical Force Measurement 


598 

Rzr 

Noise covariance 

599 

^ZR 

Bias 

600 

kzR 

Scale factor error 

601 

Rotor Roll Moment Measurement 


602 

RLr 

Noise covariance 

603 

^Lr 

Bias 

604 

•<Lr 

Scale factor error 

605 




Rotor Pitch Moment Measurement 

606 

"Mr 

Noise covariance 

607 

bMR 

Bias 

608 

"Mr 

Scale factor error 

609 




Rotor Yaw Moment Measurement 

610 

"Nr 

Noise covariance 

611 

1>Nr 

Bias 

612 

"Nr 

Scale factor error 

613 


614 

♦l 

Euler angle 

615 


Increment 

616 

®1 

Euler angle 

617 

A9 

Increment 
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Table 3.3 (Continued) 



INDEX 

SYMBOL 

DESCRIPTION 

618 

^0 ref 


619 

620 
621 

^Ic ref 
^Is ref 

Blade flapping 

Blade 1 Flapping Measurement 

622 

*^61 

Noise covariance 

623 

^61 

Bias 

624 

625 

Si 

Scale factor error 
Blade 2 Flapping Measurement 

626 

S 2 

Noise covariance 

627 

^62 

Bias 

628 

629 

'^82 

Scale factor error 
I Blade 3 FlaoDina Measurement 

630 

*’^63 

Noise covariance 

631 

^83 

Bias 

632 

633 

^83 

Scale factor error 
Blade 4 Flapping Measurement 

634 

^84 

Noise covariance 

635 

^84 

Bias 

636 

637 

•^64 

Scale factor error 
Blade 5 Flapping Measurement 

638 

*^85 

Noise covariance 

639 

^85 

Bias 

640 

641 

'^85 

Scale factor error 
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Table 3.3 (Continued) 


1 

1 


INDEX 

SYMBOL 

DESCRIPTION 



Blade 6 Flapping Measurement 

642 

^86 

Noise covariance 

643 

•^66 

Bias 

644 

645 

•^86 

Scale factor error 
Blade 7 Flapping Measurement 

646 

*^87 

Noise covariance 

647 

^87 

Bias 

648 

649 

•^87 

Scale factor error 

650 

^0 ref 


651 

^Ic ref 

Lag references 

652 

653 

^Is ref 

Blade 1 Lag Measurement 

654 


Noise covariance 

655 


Bias 

656 

657 

"cl 

Scale factor error 



Blade 2 Lag Measurement • . 

658 

R 2 

Noise covariance 

659 


Bias 

660 

661 

k^2 

Scale factor error 
Blade 3 Lag Measurement 

662 

^3 

Noise covariance 

663 

^3 

Bias 

664 

665 

'^i;3 

Scale factor error 
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Table 3.3 (Continued) 


INDEX SYMBOL 


DESCRIPTION 


Blade 4 Lag Measurement 


Noise covariance 
Bias 

Scale factor error 
Blade 5 Lag Measurement 


Noise covariance 
Bias 

Scale factor error 

Blade 6 Lag Measurement 
Noise covariance 
Bias 

Scale factor error 
Blade 7 Lag Measurement 


Noise covariance 
Bias 

Scale factor error 


Noise covariance 
Scale factor error 





Table 3.3 (Concluded) 


INDEX 

SYMBOL 

DESCRIPTION 

686 


Scale factor error 

687 


Bias 

688 

s 

Noise covariance 

689 




701 Zj Reference value for expansion variable 

702 Z 2 Reference value for expansion variable Z 2 

703 z- Reference value for expansion variable z., 

"^0 ■ 

704 z^ Reference value for expansion variable. z^ 

0 

705 Zc Reference value for expansion variable Zj. 

^0 









nwno 
129 - no 
My-iba 
Ui'J- I/O 
iby-i‘jo 
2U9-2iO 


JIU 

335 

360 

JU5-4UO 












n o o 


CH-53A NONLINEAR INPUT DESIGN RUN 100 KNOTS 
5 1 10 5 


.0 .01 
500 .01 0 

1 5 5 1 3 

• 

o 

• 

o 

.05 

• 

1 



• 

1. 

1. 

26224 .24048 .34029 

1. 

1. 

.00031554 ,00030634 

.00007458 


.3348E-064.<»698E-05 


100 . 

A 0 0 0 0 2 

I I 50 500 

y V M P 0 P PHI THETBO BOO 
MEASURE AX AY AZ QD RO P Q R P^ THET 
EXP VAR OTR 

A . ■ 


426 

C 

32.174 

427 

RHO 

.00237< 

428 

THE TAR 

0. 

4 29 

OHECA 

19.33 

430 

R 

36. 

431 

SIGHA 

0.1151 

432 

GAHHA 

10.84 

433 

CL A 

6.09 

434 

NASS 

1041. 

435 

XHUB 

0. 

436 

ZHUrt 

0. 

43 7 

IX 

41400. 

4 38 

1 Y 

206500 


Figure 3.1 Listing of Input Data 


. 000076 7957. 5 433E-05 


BC BCD BS OSD 



<•39 IZ 192500. 


^•#2 

lYZ 

0. 

A^3 

10 

4865. 


IdETA 

4865. 

<•58 

NO 

6. 

489 

AZBI AS 

32.174 

470 

RAX 

.26224 

479 

RAY 

.24048 

488 

RAZ 

. 34029 

503 

RQOaT 

.00031554 

509 

RRDOT 

.00030634 

515 

RP 

.000074580 

521 

RQ 

.000076795 

527 

RR 

. 000075433 

533 

RPHI 

.0000033346 

539 

1 

RTHETA 

CXO 

. 000004 4698 

2 

CXU 

-.000742 

3 

CXV 

♦ . 000064 9 

4 

CXW 

♦.0146 

5 

CXP 

-.000541 

6 

CXQ 

-.0428 

7 

CXR 

♦.00371 

8 

CXD FR 

-.000107 

9 

CYO 

♦ .0 

10 

CYU 

♦.00241 

11 

CYV 

-.008 7b 

12 

CYH 

♦.000334 

13 

CYP 

♦.0285 

14 

CYQ 

-.00283 

15 

CYR 

0.717 

16 

CYOFR 

♦ .00148 

17 

CZO 

♦ .0 

18 

CZU 

-.0192 

19 

CZ V 

-.0011 

2 0 

CZH 

-.115 

21 

CZP 

-.00261 


8 

16 

17 

18 

19 

20 
21 
22 
23 
11 
32 
3 3 
3 ^ 
35 


Figure 3.1 (Continued) 


-sj O' ^ uJ rv 


22 

C20 

»0.829 

36 

23 

CZR 

♦.0189 

37 

ZH 

CZDH 

.00011 ) 

38 

2*>*CL0 

.0 

66 

26 

CLU 

-.000188 

67 

27 

CLV 

-.0027 

68 

28 

CLH 

-.0015 i 

69 

29*CLP 

-.000816 

50 

30 

CLO 

.0000806 

51 

31 

CLR 

.0021 

52 

32 

CLDTR 

.000965 

53 

33 

CHO 

. .00 

61 

36 

CHU 

-.000626 

62 

35 

CHV 

.00196 

63 

36 

CHH 

•00136 

66 

37 

CMP 

-.00013 

65 

38 

CHQ 

-.1183 

66 

39 

CHR 

.000106 

67 

60 

CHOTR 

.0000096 

68 

61 

CNO 

.0 

76 

62 

CNJ 

-.00127 

7 7 

63 

CMV 

.00663 

78 

66 

CNH 

-.000159 

79 

65 

CNP 

.00265 

80 

66 

CNQ 

.00131 

81 

67 

CNR 

-.0131 

82 

68 

CNDTR 

-.0059 

83 

69 

HO 

•d 

91 

50 

HBO 

-.00568 

95 

51 

HtiOD 

-.00017) 

92 

52 

HHC 

.00533 

96 

53 

HflCD 

.0103 

93 

56 

HttS 

.000626 

9 7 

55 

HHSU 

.0021 

: 96 

56 

HIS 

.00109 

108 

57 

HTC 

.00109 

107 

58 

HTO 

-.00589 

106 


Figure 3.1 (Continued) 


YO .0 131 

60 YdO -.0022<i I 3 f, 

61 YBOD -,00lbl 132 

62 YBC .OOOTb 136 

63 Y 8 CD -.00031)2 133 

64 Y0S .005>?3 137 

65 YBSO .0039a 13 ^i 

66 YTS .00257 143 

67 YTC .0025^ 147 

68 YTO -.000157 146 

69 TO .007147 III 

70TB0 .0404 115 

71 TBOD -. 00343 112 

72 TBC .00247 H 6 

73 TBCD -.016 2 H3 

74 TBS -.0025 117 

75 T 8 SD .0193 H 4 

76 TTS -.0275 128 

77TTC -.0039 127 

78 TTO .0598 126 

79 LO .0 151 

80 LBO .000123 155 

81 LBOD .0000474 152 

82 LBC .00112 156 

83 LBCD .00224 153 

84»LBS -.00437 157 

85 L 8 Si) -.00508 154 

86 LTS -.000224 168 

87*LTC -.00031/ 167 

88 LTO .000847 166 

89 HO ,0 171 

90 MflO -.0000783 175 

91 HBOi) -.000996 172 

92 HBC .004 7 1 76 

93 HBCO .0071 173 

94 MBS .00122 I 77 

95 HUSO .00246 174 


Figure 3.1 (Continued) 


96 

HTS 

-.000531 

97 

MTC 

-.0000973 

98 

MTO 

.000598 

99 

00 

.0 

100 

080 

.0015 

lOi 

OBOO 

.000335 

102 

QBC 

.000957 

103 

08C0 

.00117 

105 

QBS 

-.0021 

105 

QBSO 

-.00175 

106 

OTS 

.000103 

107 

ore 

.00052 

108 

OTO 

-.0536 

109 

BOO 

.0 

110 

BOU 

.00059 

III 

BOV 

.00001 

112 

BOH 

. 00857 

113 

BOP 

.000275 

115 

BOO 

-.00221 

115 

BOK 

-.00221 

116 

BOBO 

-.00585 . 

117 

B0B03 

-.00735 

118 

BOBC 

.000552 

119 

BOBCO 

.000557 

120 

BOBS 

-.000585 

121 

B0BS3 

-.000325 

122 

80TS 

-.00155 

123 

BOTC 

.000283 

125 

BOOTR 

-.0000559 

125 

uoro 

.00655 

126 

BCO 

.0 

127 

BCU 

.00179 

12 8 

8CV 

-.00152 

129 

BCH 

.000751 

130 

BCP 

-.00572 

131 

BCO 

-.0105 

132 

BCR 

.00003 7 


Figure 3.1 (Continued) 



LBd 

1B7 

186 

191 

195 

192 

196 

193 

197 
195 
2 08 
207 
206 
211 
212 
213 
215 

215 

216 
217 
221 
213 
222 

219 
223 

220 
235 
233 
211 I 
2 32 

2 36 
23 7 
238 
2 39 

250 

251 

252 


133 

UC30 

.00253 

246 

134 

0C800 

.00155 

243 

135 

UCBC 

-.0013B 

247 

136 

BCBCO 

-.0103 

244 

137 

3CBS 

-.00454 

248 

138 

BCBSO 

-.0132 

245 

139 

BCTS 

-.00337 

259 

140 

UCTC 

.0747 

258 

141 

BCDTk 

-.0000153 

236 1 

142 

BCTO 

.000335 

257 

143 

BSO 

.0 

261 

144 

3SU 

-.0012) 

262 

145 

BSY 

-.00102 

263 

146 

BSH 

-.0019 

264 

147 

BSP 

-.00847 

265 

148 

BSQ 

.00674 

266 

149 

BSR 

-.000401 

267 

150 

Bsao 

-.00024 

271 

151 

BSBOD 

.00037 

268 

152 

BSBC 

.00444 

272 

153 

Bsaco 

.0127 

2 69 

154 

BSBS 

-.0014 

273 

155 

BSBSO 

-.00931 

270 

156 

'isr s 

.0047 

284 

157*BSTC 

.00317 

283 

158 

BSD IR 

.0000034 

261 1 

159 

BSTO 

-.00242 

282 

Y( 13 

LONG 

ACCEL (FT/SEC#42» 


YC 21 

KAT ACCLL (^T/SEC^^-a) 


Y( 3) 

VEKT 

ACCEL (fT/SEC»*2) 





Figure 3.1 (Concluded) 



IV. PROGRAM OUTPUT 


The sample input deck in Figure 3.1 was run on INDES. Some 
o£ the resulting output will be shown in this section. The rest 
of the output is keyed to the discussion in Section 2.2. The sum- 
mary page is shown in Figure 4.1. The namelist input variables 
are listed in Figure 4.2. The array defining one portion of the 
control input is shown in Figure 4.3. Figure 4.4 shows the begin- 
ning of the output of the optimal control sequence computed by 
INDES, while Figure 4.5 is an excerpt from the sequence at a later 
time. Figure 4.6 is a sample of the plots. which come out on the 
printer. 
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CH-S3A >IUHLINE*< INPUT DESIGN RUN 100 KNUTS 


INPUT DESIGN PROCKAH 


MRITIEN OV: 

SYSTEHS CnUTKOL. INC. 
PALO ALTUt CAL. 


THE PKUGRAN HAS SET TU THE FOLLOMING 0 1 HENS I OHS. . . 

NAS LIMIT EXCEEUE07 


lA STATES NO 

10 HEASUKEMENTS NO 

I CONTROLS NO 

72Z PARAMETERS IN MOOEL '. . NO 

SOO DATA POINTS I i TIME HISTORIES . . NO 

<i PARAMETERS TO BE lOENTIFIED 


S TIME INTERVALS 

DELTA T - .0100 

TOTAL TIME - 5.0000 SECOiOS 


3ESICH FOR CASE 1 
HITH LIMIT OH I • 5 
AND LIMIT ON J • S 

EXPANSION VARIABLES USEO IN THIS RUN... 

TIM - TAIL ROTOR 

THE FLAG ARRAYS FOR THIS RUN ARE... 

ISFLAG • I T 1 A 5 6 7 8 0 10 IT I A 9 ll 

HSFLAG • IT ) A 5 6 7 8 11 10 lA 11 15 12 

NMFLAG • 1 2 1 5 b 7 8 9 10 ll 

N/FLAG -9 


lUgure 4,1 Summary Page 


tonut 


NS 

- 

HQ 

• lOf 

HP 

- 

N3UH 

- It 

NE)(V 

• 5, 

MPfcX 

• AOO, 

NAXP 

- ^22. 

HHP 

• 5f 

N'l 

- SOOf 

NSA 

• 8^t 

NPA 

■ 60« 

NSU 

- Of 

LCV 

- 20. 

NXA 

- 0. 1, 

NrA 

- 60. 1. 

NCV 

• i5. I5t 


NR . 5, 

NIO . I, 

tCNU 


Figure 4. 2 Namelist Input 


CONTRUL NO. 1 


0. 

•tooooe -01 

0. 

0. 

0. 




0. 

0. 

0. 

0. 

0. 

u. 

Q. 

0« 

0. 

0. 

0. 

0. 

0. 

0. 

0« 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0« 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

0. 

0. 

0. 

0« 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

Q. 

0. 

0* 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

• lOOOOE-Ol 

.lOOOOE-Ol 

.lOOOOE-Ol 

•lOOOOE-Ol 

.loaooE-oi 

•lOOOOE-01 

•lOOOOE-UL 

•lOOUOE-Ol 

• lOOOOE-Ol 

. lOOOOF-Ol 

• lOOOOE-OI 

•lOOOOE-OI 

•lOOOOE-01 

.lOOOOE-01 

•lOOOOE-01 

•lOOOOE-01 

. lOOOOE’Ol 

.lOOOOE-Ol 

. lOOOOE-01 

•10000E~0l 

•lOOOOE-Ol 

• lOOOOE-Ol 

•lOOOOE-Ol 

• lOOOOE-01 

•lOOOOE-Ol 

•lOOOOC-Ol 

• lOOOOE-Ol 

•lOOOOE-Ol 

•10030E-01 

•lOOOQE-Ol 

•LQOOOE-OI 

•lOOOOE-Ol 

•10000E*01 

•10Q00E~01 

.lOOOOE-Ol 

•lOOOOE-01 

•lOOOOE-Ol 

• lOOOOE-Ol 

.lOOOOE-Ol 

• lOOOOE-Ol 

• 10000E*01 

•lOOOOE-Ol 

• iOOOOE-01 

• lOOOOE-Ol 

• lOOOOE-Ol 

• lOOOOE-Ol 

•lOOOOE'Ol 

«iOOOOE-Ol 

•lOOOOE-01 

.lOOOOE-01 

. lOOOOE'Oi 

.iOOOOE-Ol 

• lOOOOE-OL 

• lOOOOE-01 

,10090E-01 

. lOOOOE-31 

• tOOOOE-OI 

•lOOOOE-Ol 

•lOOOOE-OI 

•ioonoE-01 

* 10030E-01 

.lOOOQE-01 

•lOOOOE-01 

. tOOOOE-01 

• lOOOOE-01 

. lOOOOE-Oi 

. lOOOOE-OI 

. iOOOOE-01 

• tOOOOE-OL 

• lOOOOE-Ol 

•10000:-0l 

« lOOOOE-Ol 

. lOOOOE-Ol 

• lOOOOE-01 

• lOOOOE-Ol 

• lOOOOE-Oi 

• LOOOOE-Ol 

•lOOOOE-Ol 

•loooaE~oi 

• lOOOOE^Ol 

. lOOOOE-01 

. lOOOOE-Ol 

• iOO0OE~O| 

•IOOOOE-01 

•lOOOOE-Ol 

. lOOOOE-01 

«10000E-0I 

•lOOOOE-01 

. lOOOOE-Ol 

• lOOOOE-Ol 

• lOOOOF-Ol 

. lOOOOE-Ol 

•10030E-0L 

• lOOOOE-01 

•lOOOOE-Oi 

•IOOOOE-01 

• lOOOOE-0 1 

• lOOOOE'Ol 

• lOOOOE-Ot 

0. 

0. 

0. 

0. 

0 . 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0, 

0. 

0 . 

0, ' 

0. 

0. 

0. 

0. 

0« 

0. 

0. 

0. 

0, 

0. 

0. 

0. 

0. 

0. 

0. 

u« 

0. 

0. 

0. 

0. 

0. 

0. 

0« 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 



Figure 

4.3 Control Input 









Figure 4.3 (Continued) 


SYSTEM SIHULATION iflTH OPTIHAL INPUTS 


TIHEISECI CUNTRILS — 


KEASUREHENTS — 


0.000 0 . 

•010 0« 

.020 0. 

•030 0. 

00 • 0 0 • 
• 0*>0 0 . 

•000 0. 

.0^0 0. 

•oao 0 . 


f.67^llE-02 0, 0. 
7.67411E-02 0. 0. 
7.67AIIE-02 0. 0. 
7.67SI1E-02 0. 0. 
7.6741IE-02 0, 0, 
7.67SIIE-02 0. 0. 
7,67^lIE-02 0, 0. 
7.67AilE>02 0. 0. 
7.67^116-02 0. 0. 


0. 

0. 

-^•169006-02 

0. 

0 . 

0. 

-1.19107E-OI- 

9.63108E-06 

0. 

0 . 

-1.5A500E-01- 

7.66636C-0S 

0. 

0. 

-1.49672E-01- 

2.53575E-0S 

0. 

0. 

-1.0858^6-01- 

5.85016E-05 

0. 

0. 

-3.718396-02 
1. 106656-03 


0 . 

1.92962: tOl 

0 . 

0. 

1.55337i»Ol 

7.I5520E-O5 

0. 

1.00550=»01 

5.39015E-03 

0. 

6.30126E»00 

5.95292£-03 

0. 

3«25555:»00 
1. 16765E-03 

0. 

9.51579E-01 

1.32657E-02 


0 . 0 . 

0. -6.60620H-01 

5.755556-02 3.985916-02 
1.855626-03 
0 . 0 . 

0. -1.606556*00 

I.60125E-OI 8.016156-02 
2.815866-03 
0 . 0 . 

0. -1.976056*00 

2.67I67E-01 1.359276-01 
5.028566-03 


0. 0. 0. 

9.12053E-01 1.3 126 IE *00-6.100396-03 

0 . 0 . 0 . 

0. 0. 0. 

9.953916-01 7.180856*00 5.70937E-01 
2.9O21OE-O3-0.39559E-O5 1.01765E-06 

0 . 0 . 0 . 

1.69 7006 *00 1.096 33E *0 I 1. 10 55 06 *00 
1. 1 302 5E -02-2. 25 3 I IE -0 3-2. 5 3 756E -05 

0 . 0 . 0 . 

2.915 7 76 *00 1.315576*01 1.55 295E*00 
2.5 75626 -02-3. 79655E-03-7. 9737 16-05 

0 . 0 . 0 . 

5.528596*00 1.522356*01 1.909756*00 
5.2 1 3826 -02-5. 11 576E -03- 1 . 22 2 726 -05 

0 . 0 . 0 . 

6.512256*00 1.559316*01 2.17137E«00 
6.262 32E -02-5. 06 707E -03-7. 1 161 2E -05 

0 . 0 . 0 . 

8.552296*00 1.56295E*01 2.33955E*00 
8. 525 356 -02 -5. 781 62E -0 3 1.8 35 76F-05 

0 . 0 . 0 . 

l.0505^E*01 1.563516*01 2.520006*00 
1.09125E-01-5.66336E-03 7. 716026-05 

0 . 0 • 0 . 

1.250356*01 1.58519E*0l 2.52557C*00 
1. 3 35 15E-01-2.39 100E-0 3 1.815 76E -0 3 


Figure 4.4 Beginning of Optimal Control Sequence 


0 


0 


VO 


.090 

• too 

• no 

• 120 
• MO 


0. 


0. 


0. 


0. 


• MO 0. 


7.67«tlE-02 0. 
7.67^11E-02 0* 
7.67^LIE-02 0. 
7.67^116-02 O. 
7.67^llE-02 0. 
/.#,/<•! i« -o > n 


0. 


0. 


0. 


0. 




0 . 0 . 0 . 0 * 0 . 

0, -l.a66l6E*00 l«^3565E*0l 1.542236*01 2.365746*00 
4.070406-01 2.04l96i-0l 1.574196-01 1.00529E-03 3.518436-03 
5.401246-03 

0. b. 0. 0. 0. 

0. -1,376696*00 1.600476*01 1. 640426*01 2.255286*00 
5.242316-01 2.074115-01 1.005616-01 5.751396-01 5.965476-03 
7.171006-03 

0 . 0 . 0 . 0 . 0 . 

0. -6.130135-01 1.740736*01 1.600296*01 2.10/066*00 
6.315576-01 3.035495-01 2.024026-01 1.15395E-02 9.310046-03 
9.009206-03 

0 . 0 . 0 . 0 . 0 . 

0. 1.257265-01 1.054106*01 2.022806*01 1.916696*00 
7.262646-01 4.912466-01 2.226396-01 1.014036-02 1.167766-02 
1.121406-02 

0 . 0 . 0 . 0 . 0 . 

0. 1.351725*00 1.940206*01 2.209236*01 1.751/66*00 
0.060006-01 6.009206-01 2.410976-01 2.601406-02 1. 917406-92 
1. 351026-02 


Figure 4.4 (Continued) 


Cn 

O 


1.960 

1.9/0 

1.900 

1.990 

^•000 

^.0l0 

2.020 

2.030 

2 * 0^0 

2.050 


0. 5.99012E-02 0. 0. 


0. 5.99032E-02 0. 0. 

0. >l.65025E*03 0. 0. 


0. -1.65025E-03 0. 0. 


0. -1.65025E-03 l.92962E»0l 9.I20^3E-01 

0. 0. i.4^337E«0l 9.95391E-01 

1.01/65E-06 9.63108E-06 1.005^0E»0i 1.69708E*00 
‘-2.%J/56E-05 7.66638E-05 6.30126E«00 2.9l^77E«00 
-7.37J71E-05 2.53575E-0V 3.235S5E«00 ^.52839E»00 
-l.22272E-0^ 3;850lBE-03 9.^1579E-01 6.41225E»00 


-“I.tiicioc-ut u.j:>iJlt»OU 1 .20 H 2 .65 8 2,91 -*UL 

-1.14385E-0I 

0. 0. 0. 0. 

0. -7.1972 H -01 1 .33961E ♦02-2.3^ 7 62 E »01 

-3.99083C-02 8.36893t*00 I . 29 600E -01-6.66 6 79E -0 I 
-I.IJ006E-01 

0. 0. 0. 0. 

0. -5.33593E-01 i . 6 3869E ^02-2.693 18E *0 1 

-6.91715E-02 6.36626EtOO 1 . 30698E-0 1-6. 67026E -0 1 
-1.11229E-01 


l.23662L»Ul 

0. 

1.03029E-01 

1.26282E»01 

0 . 

1.133696-01 

1.271226*01 


0 . 

0 . 

3.36076E-02 

-1.09071E-01 

0 . 

0 . 

I.2i/10E-01 

-1.06370E-01 

t.3126lE»00' 

0. 

2.03319E-01 
-L.03766E-0I 
7.10083E»00 
-8.39569E-06 
2.80930E-0I 
-1.00770E-01 
l.09633E*01 
-2.263116-03 
3.676986-01 
-9.75760E-02 
1.313676*01 
-3. 796356-03 
6.066666-01 
-9.623696-02 
1.622336*01 
-3.115766-03 
6.307786-01 
-9.078176-02 
1.639316*01 
-5.867076-03 
6.836126-01 
-8.723566-02 


0» 0« 0* 0* 

l.aa362E*00 t.43036E»02-).l6S19E*01-1.0 )<)I5E-0I 
a.33939£*00 t . 309^aE -01-^.671 7ir-0 t l.7796IE«01 

0. 0. 0. 0. 

%.729I25»00 1.42S91E*02-3.63323E*01-S.3%9‘I2E-01 

а. 336&5£»00 t.27707E-0t-^.66392E-0l 1.267991 »01 

б. l0039£-03-^.la990E-02 0. 0. 

9. 0 79955 *00 1 . 4 15a^E *02-3.91 70aE *01-9.09 32 3C-01 
a.33417E*00 1.20525E-01-4.64751E-01 1.2963/E*01 

5. 7a937E>01-1.19107E-01-7.35928E-05 2.90210E-03 
l.la033E»01 1. 40095E *02-5. O^a 19E *01-1. 20442E *00 
B.32a29:*00 1.097a2E-01-5.b2310E-01 1.30575E*01 

1. 10550' *00- 1.55 5a0E -01-5. 3901 5E -03 1.1 3B25E -02 
1.3a712E*01 1.38163E*02-5.0b713E*01-1.5159 3E *00 
a.31592£*00 9.bbl96E-02-5.59160E-01 1.31310E*01 

1. 552955 * 00- 1. 595 72E -01-5. 95292E -0 3 2.5 75ft2E-02 
1.53 32b: *01 1 . 35aa9E *02-5.01053E *0 1-1. 55 /22E *00 
8.29582E*00 a.l755bE-02 5.553B9E-01 1.32151E*01 

1. 909755 *0n-1.0a5B5E-0 1-1. 187b5E -03 5.213B2E-02 
1.62 150: *01 1.33 35 6E *02-3.703 155 *01 -1.6050 7( *00 
8.263556*00 6. 59251E -02-5.5 1 103E -01 1.32975E*01 

2.17137£*00-3.71B39E-02 1.32657E-02 6.26232E-02 
1. 656755*01 1. 30656E *02-3. 7659 75 *01-1.595606 *00 
0.221256*00 5.90735E-02-5.56512E-01 1.31001C*01 


Figure 4.5 Later in the Optimal Control Sequence 





2 * 0/0 


2*oao 


2.000 


2.100 


2.110 


-/.ll!il2E-05 l.l065^E-03-6.60620E-0l 8.^^220£»00 
1.6iS/6E-0^ l.8^^62E-03-1.506^5Et00 l.0505^E*0l 
/./lb02£-0^ 2 . 8l5a6E-01-l.9/605E »00 l.2503AE»0l 
l.81^76E-01 ^•02856E-03-l.865lf»E^00 l.^3^65E*0l 
3.5l8^1E-0) ^ .^832^E-01-I.17669f>00 l.600^7E*0l 
5.00aS7E-0 3 7.1 738 flb -0 3-6. I 38 1 IE -01 1. 7^07 3E *01 


l•^629^£♦01 
-5.78162E-03 
5.08327E-01 
-8.36110E-02 
l.^634lE*0l 
-^.66336E-03 
5.19719E-01 
-7.99171E-02 
1.^8Sl9E*0l 
-2.39l0n£-0i 
5.20158E-01 
-7.61599E-02 
1.5^223E»01 
1.08529E-03 
^•1056SE-01 
-7.23^366-02 
1.648626*01 
5.75139E-03 
4.92200E-01 
-6.847291-02 
1.808296*01 
1. 1939SE-02 
4.664 3lir-0l 


2.339996*00 9.79449E-02 3.9899IE-02 8.929396-02 
1.64924: *01 1.2 7840E *02-3.6 1 122E *01- 1. 9292 9E *00 
8. 16 773c *00 3.42291E-02-4.414346-01 1. 346236*01 

2.420886*00 1.681246-01 6.016196-02 1.091246-01 
1.994 106*01 1.250136*02-3.492 776*01-1.406966 *00 
a.l0329:*00 1.993196 -02-4. 16284E -01 1.39438C*Ol 

2.42947:*00 2.H7167E-01 1.14927E-01 1. 33419E-01 
1.910696*01 1.22 2 3 36*02-3.296946 *01 -1.290306*00 
8.028696*00 6.217986-01-4.310766-01 1.362476*01 

2.369745*00 4. 078406-01 2.041966-01 1.974196-01 

1.402295 *01 1.1999 76 *02-3.146226 *01-1.069746 *00 
7.3 44 846*00-9.3819 76-03-4.299146-01 1.3 704 76*01 

2.299285*00 9. 242316-01 2.874116-01 1.809616-01 

1 .2 79 79c *0 I 1.1 70 346 *02-3.002 716 *0 1-a. 63 16 36-01 
7.8 93196*00-1.9017 36-02-4.208946-01 1.3 78 386 *01 

2.107866*00 6.319976-01 3.839496-01 2.024026-01 

1. I 3 72 8; *01 1.14 7006 • 0? -2 . «6 49 U ♦ 0 1 -6. 9 1 4846 -0 1 
7. 79 90 If • 00-2.26 I 4H1 -02-4 . I6099f -0| I . 34M *ir • 


Figure 4.5 (Continued) 


\ 


A LONG A -3.00E^0l -2.00S^0l -l.OOE^Ol 0. L.OOE^Ol 2.00E*0l J.OOE^Ol 

TIHE(SEC) 

0.00000 • : A 

• 0 LO 00 . * A 

•02000 • : A 

•03000 • : A 

•OAOOO • : A 

•09000 • :A 

•06000 • a: 

•07000 • A : 

•08000* A* 

•OOOOO • A ! 

•lOOOO • a: 

• UOOO • a: 

•12000 • A 

•13000 • :a 

•lAOOO • : A 

•15000 • : A 

•16000 • : A 

.17000 • : A 

•18000« *A 

•17000 • : A , 

•20000 • £ A 

•21000 • : A 

•22000 • : A 

•23000 • : A 

•2^000 . : A 

•25000 . : A 

•26000 • : A 

.27000 . * : A 

•28000 • : A 

•27000 • : A 

•30000 • : A 

•31000 • : A 

•32000 • : A 

•33000 • : A 

•36000 • : A 

•35000 • i A 

•36000 • : A . 

•37000 • : A 

•38000 • : A 

•37000 • : A 

•60000 • : A 

•61000 • ; A 

•62Q00 • : A 

•6)000 • : A 

•66000 . : A 

•65000 • :A 

•66000 • :A 

•67000 . :A 

•68000 • :A 

•67000 • £A 

•50000. £A 

.51000 . :A 

.52000 . A 

.53000 • A 

•56000 . A 


Figure 4.6 Sample of Printer Plots 
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